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FL0WF1ELD  OF  A  LIFTING  HOVERING  ROTOR— A  NAVIER-STOKKS  SIMULATION 


SUMMARY 


The  viscous,  three-dimensional  flow  field  of  a  lifting  helicopter  rotor  in  hover  is  calculated 
by  using  an  upwind,  implicit,  finite-difference  numerical  method  for  solving  the  thin  layer  Navicr- 
Stokes  equations.  The  induced  effects  of  the  wake,  including  the  interaction  of  tip  vortices  with 
successive  blades,  are  calculated  as  a  part  of  the  overall  flowfield  solution  without  using  any  ad  hoc 
wake  models.  Comparison  of  the  numerical  results  for  the  subsonic  and  transonic  conditions  show 
good  agreement  with  the  experimental  data  and  with  lire  previously  published  Navier-Stokes  calcu¬ 
lations  using  a  simple  wake  model.  Some  comparisons  with  Euler  calculations  are  also  presented, 
along  with  some  discussions  of  the  grid  refinement  studies. 

INTRODUCTION 


The  accurate  numerical  simulation  of  the  flowfield  of  a  lifting  helicopter  rotor  continues  to 
be  one  of  the  most  complex  and  challenging  problems  of  applied  aerodynamics.  This  is  true  in 
spite  of  the  availability  of  the  present  day  supercomputers  of  Cray-2  class  and  improved  numerical 
algorithms.  However,  many  advances  have  been  made  to  date  with  the  use  of  simpler  set  of  equa¬ 
tions  of  fluid  motion,  such  as  the  potential  flow  equations,  to  model  these  complex  flowficlds.  The 
equations  have  been  simplified  by  coupling  the  solution  scheme  with  an  empirical  wake  model  to 
bring  in  the  influence  of  the  vortex  wake.  Solution  schemes  that  use  this  idea  are  often  grouped 
under  methods  using  wake  models  and  encompass  the  potential  flow  (refs.  1-5),  the  Euler  (refs.  6- 
8)  and  the  Navier-Stokes  methods  (refs.  9-12).  In  contrast  to  these  methods  that  use  ad  hoc  wake 
mode1**  there  arc  methods  that  compute  the  essential  details  of  the  induced  effects  of  vortex  wake 
as  a  f  the  overall  flowfield  solution.  These  are  called  the  wake  capturing  schemes  and  have 
been  nstrated  for  the  solutions  of  the  potential  flow  (ref.  13),  the  Euler  (refs.  14-17)  and  the 
Navier-Stokes  equations  (ref.  18). 

The  basic  assumptions  of  potential  flow  methods  restrict  their  application  to  low  supercritical 
speeds  without  the  use  of  entropy  corrections.  Despite  this  feature,  the  potential  flow  methods,  cou¬ 
pled  with  a  wake  modeling,  have  been  very  useful  in  the  industrial  environment  for  design  analysis 
(refs.  2-5).  On  the  other  hand,  the  Euler  equations  contain  the  essential  physics  to  describe  convec¬ 
tion  of  vortical  flows  and  do  not  have  the  restriction  on  the  Mach  number.  But  they  still  lack  the 
essential  ingredients  to  model  the  separated  flows  and  inviscid-viscous  interactions  associated  with 
shock-induced  separated  flows.  Nevertheless,  the  Euler  methods  have  been  used  to  model  these 
complex  vortical  flows  by  coupling  with  wake  models  (refs.  6-8)  as  in  potential  flow  methods.  But 
the  major  drawback  of  these  methods  is  that  they  have  proven  to  be  more  expensive  in  comparison 
to  the  potential  flow  methods. 

Even  Navier-Stokes  methods  (refs.  9-12)  have  been  coupled  with  wake  models  to  calculate 
these  complex  flows.  Although  these  methods  capture  viscous  effects  adequately,  they  remain 
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limited  by  the  wake  modeling,  which  tends  to  be  restricted  to  simple  geometries  and  planfonns.  In 
general,  a  major  disadvantage  of  these  methods  which  use  wake  modeling  is  that  the  technique  of 
prescribing  u.  wake  has  to  be  specialized  for  each  blade  shape  and  planform  and  therefore  cannot 
easily  handle  arbitrary  blade  shapes  with  twist  or  taper. 

Therefore,  the  weak  link  in  the  above  wake-coupled  methodologies  has  been  the  wake  mod¬ 
eling.  In  contrast  to  the  methods  using  wake  modeling,  several  schemes  have  attempted  to  capture 
the  wake  and  its  effect  as  a  part  of  the  overall  solution  scheme.  lltcse  methods  range  in  complexity 
from  potential  flow  methods  (ref.  13)  to  a  Navicr-Stokcs  method  (ref.  IB).  All  of  these  inviscid 
methods  (refs.  13-17)  utilize  finite-volume  formulation  for  the  solution  method.  Of  these  different 
wake  capturing  schemes,  the  potential  flow  scheme  of  Ramaehandran  ct  al.  (ref.  13)  appears  to  be 
tne  most  accurate.  All  of  the  Euler  methods  appear  to  compute  the  flow  in  the  tip  region  reason¬ 
ably  well.  However,  the  inviscid  methods  still  lack  the  ability  to  capture  accurately  the  formation 
of  a  tightly-braided  tip  vortex  structure,  and  therefore,  the  accuracy  of  the  computed  wake  and 
tip-vortex  core  may  be  questionable. 

The  purpose  of  this  study  is  to  develop  a  calculation  method  for  the  solution  of  Navicr-Stokes 
equations  for  the  complete  flowfield  of  a  lifting  rotor,  including  the  wake  and  its  induced  effects. 
'Hie  vortex  wake  and  its  effects  arc  captured  as  a  part  of  the  complete  llowlield,  and  thus  no  arbitrary 
inputs  are  necessary  to  describe  the  wake.  Although  this  is  not  very  different  in  concept  from 
the  Euler  wake-capture  schemes  discussed  above,  the  Navicr-Stokcs  approach  was  needed  for  the 
following  reasons:  1)  better  tip-flow  simulation,  which  involves  resolving  the  blade-lip  separation 
and  the  formation  of  a  concentrated  tip  vortex,  2)  accurate  simulation  of  strong  viscous-inviscid 
interaction  involving  shock  induced  separation  at  high  blade  tip  speeds  and  high  collective  pitch 
conditions,  and  3)  future  modeling  of  retreating  blade  and  dynamic  stall  regimes  in  forward  flight. 

The  numerical  code  used  in  this  study  is  an  improvement  of  the  version  that  was  developed 
previously  in  related  studies  with  wake  modeling  (ref.  9).  One  fundamental  difference  of  the  new 
numerical  scheme  is  the  use  of  Roe’s  upwinding  in  all  three  directions  (ref.  19).  litis  feature,  cou¬ 
pled  with  a  simplified  left-hand-side,  has  produced  an  efficient  and  accurate  numerical  scheme. 
These  additional  changes  in  the  Navicr-Stokcs  algorithm  are  based  on  some  of  the  numerical  pro¬ 
cedures  described  in  reference  20  and  will  be  described  briefly  in  the  following  sections. 

The  first  author  (GRS)  would  like  to  acknowledge  the  support  of  this  research  by  the  U.S. 
Army  Research  Office  under  Contracts  DAAL03-88-C-0006  and  DAAL03-90-C-0013.  Computa¬ 
tional  time  was  provided  by  the  Applied  Computational  Fluids  Branch  of  NASA  Antes  Research 
Center. 


GOVERNING  EQUATIONS 


The  governing  differential  equations  are  the  thin  layer  Navicr-Stokes  equations.  These  can  be 
written  io  conservation-law  form  in  a  generalized  body-conforming  curvilinear  coordinate  system 
as  follows  (ref.  21): 

drQ  +  diE+dnF+d(id=’l-d<S  (1) 
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where  t  *  t,  (  =  ((x.y.z,*)*  t?  *  tj (x,y,ztt)  and  {  *  £(x,y,z,0-  The  coordinate  system 
(x,y,z,0  is  attached  KHhe  blade  (see  fig.  1).  The  vector  of  conserved  quantities  Q  and  the  inviscid 
flux  vectors  B ,  F,  and  5  are  given  by 
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where  H  *  (e  +  p)  and  c  *  0  or  1  for  the  Euler  or  the  Navier-Stokes  equations,  respectively.  In 
these  equations,  U,  V,  and  W  are  the  contravariant  velocity  components  defined,  for  example,  as 
U  =  &  +  (*u  +  £yv  +  £tw.  The  Cartesian  velocity  components  are  represented  by  u,  v,  and  w 
and  the  density,  pressure,  and  the  total  energy  per  unit  volume  by  p,  p,  and  e,  respectively.  The 
characteristic  length  and  velocity  scales  are  the  rotor  blade  chord  and  the  ambient  sound  speed,  and 
p  and  p  are  nondimensionalized  by  their  respective  ambient  values.  The  quantities  £z»  £y.  6*  (t< 
etc.  are  the  coordinate  transformation  metrics  and  J  is  the  Jacobian  of  the  transformation.  For  the 
thin  layer  approximation  used  here,  the  viscous  flux  vector  5  is  given  by 
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rot  -  cj +  cj +  c? 

m2  =  C*u< +  <vy<  + 

raj  =  i(uJ  +  #l  +  (o2)< 

where  Re  is  the  Reynolds  number,  Pr  is  the  Prandtl  number,  7  is  the  ratio  of  specific  heats,  and  a 
is  the  speed  of  sound.  The  fluid  pressure,  p  is  related  to  the  conservative  flow  variables,  Q,  through 
the  nondimensional  equation  of  state  for  a  perfect  gas, 

P=  (7- !){e- ^(u2  +  u2 +  ui2)  }  (4) 


For  turbulent  viscous  flows,  the  viscosity  coefficient  p  in  S  is  computed  as  the  sum  of  pj  +  pt 
where  the  laminar  viscosity,  pi,  is  estimated  using  Sutherland’s  law  and  the  turbulent  viscosity,  p* , 
is  evaluated  using  the  Baldwin-Lomax  algebraic  eddy  viscosity  model  (ref.  22). 
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NUMERICAL  ALGORITHM 


A  finite-difference,  upwind,  numerical  algorithm  is  developed  for  the  helicopter  rotor  appli¬ 
cations.  The  evaluation  of  the  inviscid  fluxes  is  based  on  an  upwind-biased  flux-difference  split¬ 
ting  scheme  for  the  right-hand  side  while  an  LU-SGS  (Lower-Upper  -  Symmetric  Gauss-Scidel) 
scheme,  suggested  by  Jameson  and  Yoon  (refs.  23  and  24),  is  used  for  the  implicit  operator.  The  van 
Leer  MUSCL  (monotone  upstream-centered  scheme  for  the  conservative  laws)  approach  (ref.  25) 
is  used  to  evaluate  the  conservative  variables  to  obtain  the  second-  or  third-order  accuracy  with 
flux  limiters  so  as  to  be  TVD  (total  variation  diminishing).  The  upwind-biased  scheme  used  on  the 
right-hand  side  was  originally  suggested  by  Roe  (ref.  26 )  and  later  extended  to  three-dimensional 
flows  by  Vatsa  ct  al.  (ref.  19).  The  chief  advantage  of  using  upwinding  is  that  it  eliminates  the  ad¬ 
dition  of  explicit  numerical  dissipation  and  is  known  to  produce  less  dissipative  solution  (ref.  19). 
This  feature,  coupled  with  a  fine  grid  description  in  the  tip  region,  increases  the  accuracy  of  the 
wake  simulation.  A  similar  algorithm  was  used  in  the  finite-volume  Euler  scheme  of  reference  17 
to  investigate  the  exact  same  problem  studied  here. 


The  space-discretizcd  form  of  the  differential  equation,  equation  1,  is 

a*  M i h-i 

dtQ~~  ^ 

1  -  Si-i 

AC  *  Rc  AC 


(5) 


where  ; ,  Jfc,  and  l  correspond  to  the  £,  tj,  and  C  coordinate  directions,  respectively. 

The  application  of  Roc’s  upwinding  (ref.  19)  to  the  numerical  flux  of  the  inviscid  terms  results 
in  the  locally  one-dimensional  form  and  can  be  written,  c.g.,  in  the  £  direction,  as 


-  -  Q,.)] 


where  A  is  the  Roc-averaged  Jacobian  matrix  and  Qi,  and  Qn  are  the  left  and  right  state  variables. 
The  scheme  degenerates  to  the  first-order  accuracy  if  Qi  =  Qj  and  Qji  =  Q/+i .  Higher-order 
schemes  can  be  constructed  from  a  one-parameter  family  of  interpolations  for  the  primitive  vari¬ 
ables,  p,  u,  v,  u»,  and  p.  For  example, 

PI  =  {1  +  ^[(1  -iOV  +  O  +  kJADp/ 

Pr  =  {|-%Ll<1  +  <°V+(l-/t)A]}pjtl 


where  V  and  A  are  backward  and  forward  difference  operators,  and  k  is  a  parameter  that  controls 
the  construction  of  higher-order  differencing  schemes.  For  example,  to  construct  the  third-order 
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scheme  in  the  present  method,  k  = 
V>  is  calculated  as 

V»/  = 


y,  Koren’s  differentiable  limiter  (ref.  27)  is  used.  The  limiter 
3VpyAP/-f  c 

2(  A  pj  —  Vpj)1  +  3  Vp/Ap;  +  c 


where  a  small  constant,  typically  c  *  10"*,  is  added  to  prevent  the  division  by  zero.  Similar 
formulae  are  used  for  the  other  primitive  variables.  The  viscous  flux  terms  are  discretized  using 
second-order  central-differencing  (ref.  21). 

The  time  marching  integration  procedure  uses  the  LU-SGS  method.  The  details  of  this  scheme 
are  described  elsewhere  (ref.  20).  The  final  form  of  this  algorithm  can  be  written  for  a  first-order 
time  accurate  scheme  as 

LDUAQ*  *  -AtRHSn  (8) 


where 
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D  =  t/  +  &t(A*  -A~  +  B-B~  +  C~  CHI;, 

U  =  I  +  4ii*|w  +  &t&(A-  +  Atg‘|M,i  +  AlA,g-  +  A  (0*1;, tli  +  A(A<C" 

where  At  is  the  time  step,  RHS  represents  the  discretized  steady  state  terms,  e.  g.,  equation  5,  and 
nrefers  to  the  current  time-level.  Also,  A*  -  y(A  +  <T(),  A~  =  y(/l  -  aj).  cr^  =  |C7|  +  ar(  +  e, 

e  =  0 .01  typically,  and  +  (v2  +  £*2.  As  a  result  of  the  simplified  form  of  the  Jacobian 

terms,  e.g.,  A+,  the  block  diagonal  matrix  D  reduces  to  a  scalar  diagonal  matrix.  Thus  this  method 
requires  only  two  (one  forward  and  one  backward)  sweeps  with  scalar  inversions  and  leads  to  less 
factorization  error.  Lastly,  additional  source  terms  have  been  introduced  to  account  for  the  rotation 
of  the  blades  because  of  the  blade-fixed  coordinate  system  used  here. 

The  present  numerical  scheme  uses  a  finite-volume  method  for  calculating  the  metrics.  The 
chief  advantage  of  such  a  formulation  is  that  the  metrics,  including  the  time  metrics,  can  be  formed 
accurately  (ref.  28),  and  this  approach  captures  the  free-stream  accurately  (ref.  17).  To  be  compati¬ 
ble  with  the  present  finite-difference  numerical  scheme,  the  metrics  are  evaluated  at  the  grid  nodes 
instead  of  the  cell  centers  of  a  standard  finite-volume  method  (ref.  20).  Also,  the  time  metrics  arc 
evaluated  in  the  same  manner  as  in  a  finite-difference  scheme,  which  is  less  expensive  compu¬ 
tationally  than  rigorous  evaluation  of  the  time  metrics.  However,  free-stream  subtraction  is  then 
required  to  restore  accuracy  to  the  time-metric  terms. 

The  flowfield  of  a  hovering  rotor  is  initially  quiescent  (ref.  29)  and  the  evolution  of  the  flow- 
field  is  monitored  as  the  blade  is  set  in  motion.  To  take  advantage  of  the  quasi-steady  nature  of 
the  hovering  rotor  flowfield,  a  locally-varying  time  step  is  used  in  the  integration  procedure  to 
accelerate  convergence,  as  suggested  in  reference  30. 
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Grids  and  Boundary  Conditions 


Body -con forming  grids  were  generated  for  the  rotor  blades  using  an  elliptic  solver.  Because  of 
the  cylindrical  nature  of  the  flow  of  a  hovering  rotor,  a  C-H  cylindrical  grid  topology  was  chosen, 
as  in  reference  17.  In  contrast  to  the  experimental  model  rotor  that  has  a  square  tip,  the  present 
numerical  scheme  approximates  this  as  a  bevel  tip  because  of  the  M-topology  of  the  grid  in  that 
direction  (sec  ref.  31). 

The  standard  viscous  grids  used  here  had  217  grid  points  in  the  wrap  around  (along  the  chord) 
direction  with  144  points  on  the  body,  71  points  in  the  spanwisc  (radial)  direction  with  55  points 
on  the  blade  surface,  and  61  points  in  the  normal  direction.  The  grid  was  clustered  near  the  leading 
and  trailing  edges  and  near  the  tip  region  to  resolve  the  tip  vortex.  It  was  also  clustered  in  the 
normal  direction  to  resolve  the  viscous  flow  near  the  blade  surface.  There  are  about  15  points  in 
the  boundary  layer  with  a  spacing  of  the  first  grid  point  from  the  surface  equal  to  5x  10"5  chord 
(that  translates  to  a  y+  =  0(1)).  A  coarse  grid  was  constructed  from  this  fine  viscous  grid  by 
removing  every  other  point  in  all  three  directions.  The  inboard  plane  near  the  axis  of  rotation  was 
located  at  a  radial  station  equal  to  one  chord.  The  grid  outer  boundaries  were  set  at  8  chords  in  all 
directions.  The  same  grids  were  used  for  the  Euler  calculations. 

Figure  1  shows  the  coarse  grid  that  was  used  in  the  computations.  Because  of  the  symmetry 
of  the  hovering  flow  and  the  periodic  boundary  condition  described  below,  the  calculations  could 
be  performed  for  only  one  blade.  Figure  la  shows  the  cylindrical  nature  of  the  grid  in  the  plane 
of  the  rotor,  and  figure  lb  shows  an  isometric  view  of  the  grid  boundary  for  a  single  blade.  For 
clarity,  the  figure  shows  only  the  blade  and  side  boundaries.  The  bottom  surface  and  other  grid  line 
are  omitted.  Also  shown  is  the  coordinate  system,  where  x  is  in  the  chordwise  direction,  y  is  in  the 
radial  direction,  and  z  is  in  the  normal  direction.  The  blade  motion  is  counterclockwise. 

All  the  boundary  conditions  are  applied  explicitly.  The  radial  inboard  and  far-ficld  boundaries, 
as  well  as  the  upper  boundary  of  the  cylindrical  mesh,  are  updated  by  means  of  a  characteristic- 
type  boundary  condition  procedure,  although  the  Roe's  upwinding  used  in  the  numerical  procedure 
would  otherwise  treat  the  boundaries  in  a  1-D  characteristic  sense  anyway.  At  the  wall  a  no-slip 
boundary  condition  is  used  for  the  viscous  calculations.  The  Euler  calculations  use  an  extrapolation 
of  the  contravariant  velocities  at  the  surface.  The  density  at  the  wall  is  determined  by  a  zeroth- 
order  extrapolation.  The  pressure  along  the  body  surface  is  calculated  from  the  normal  momentum 
relation  (sec,  for  example,  ref.  21).  The  total  energy  is  then  determined  from  the  equation  of  state. 

To  capture  the  information  in  the  wake  region  of  the  blade,  a  periodicity  condition  is  used 
to  swap  the  information,  after  interpolation,  at  the  front  and  back  boundaries  of  the  cylindrical 
grid  topology  (see  fig.  lb).  This  is  also  done  in  an  explicit  manner.  At  the  bottom  boundary,  the 
scene  of  the  far-field  wake,  an  approximate  condition  based  on  the  normal  velocity  is  used.  For  an 
outflow  condition,  all  conserved  flow  quantities  are  extrapolated  from  the  grid  interior  except  for 
the  energy,  which  is  calculated  by  prescribing  the  free-stream  pressure.  For  inflow  at  this  boundary, 
the  free-stream  (ambient)  conditions  are  specified. 
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RESULTS  AND  DISCUSSION 


The  test  cases  considered  in  this  study  correspond  to  the  experimental  model  Stover  test  con¬ 
ditions  of  Caradonna  and  Tung  (Ref.  29).  The  experimental  model  consists  of  a  two-bladed  rigid 
rotor  with  rectangular-planform  blades  with  no  twist  or  taper.  The  blades  are  made  of  NACA  0012 
airfoil  sections  with  an  aspect  ratio  of  6.  Three  experimental  conditions  were  chosen  from  among 
the  dau:  1)  tip  Mach  number  A/t*>  =  0.44,  collective  pitch  9  *  8*,  and  the  Reynolds  number  based 
on  the  tip  speed,  Re*  1.92x  106;  2)  Af,,p  *0.877,  9  *  8*  and  Re*  3.93x  10*;  and  3)  Affip  *0.794, 
0  =  12*  and  Re  =  3.55x10*. 


Fine  Grid  Navier-Stokes  Results 

Surface  pressures  are  shown  in  figures  2-4  for  the  three  experimental  conditions  considered. 
These  calculations  were  done  on  a  fine  grid  consisting  of  nearly  one  million  points.  Figure  2  shows 
the  surface  pressures  for  the  conditions  of  Af„p  *  0.44, 6  =  8°,  and  Re  =  1 .92x  10* .  In  this  figure,  the 
present  calculations  are  compared  with  the  experimental  data  of  reference  29  and  the  results  from  a 
previous  Navier-Stokes  calculation  that  used  a  simple  wake  model  (ref.  9).  The  present  calculations 
agree  well  with  the  experimental  data  for  all  radial  stations.  There  are  some  improvements  in  the 
results  at  y/R  =  0.50  and  0.96  over  the  previous  results  from  reference  9.  It  should  be  pointed 
out  that  the  calculations  of  reference  9  used  a  0-0  grid  topology  with  nearly  700,000  grid  points 
having  a  grid  clustering  similar  to  the  present  grid. 

Figure  3  shows  a  comparison  of  surface  pressures  for  the  condition  of  A/t,>  *  0.877, 9  *  8s  and 
Re  =  3.93  x  1 0* .  At  this  transonic  flow  condition,  the  present  calculations  show  excellent  agreement 
with  the  experimenni  data  for  all  radial  stations.  In  contrast  to  the  calculations  of  reference  9,  the 
present  results  show  shock  locations  and  shapes  that  are  well  captured.  The  inboard  regions  of  the 
flow  are  also  predicted  mere  accurately;  this  indicates  that  the  present  computed  wake  is  superior 
to  the  approximate  wake  model  of  reference  9. 

Figure4  shows  a  comparison  of  surface  pressures  for  the  condition  of  A/„-p  =0.794, 9  =  12°  and 
Re  -  3.55  x  10* .  Because  of  the  high  collective  pitch,  this  case  is  more  severe  in  terms  of  the  shock 
strength  and  shock-induced  boundary  layer  separation,  even  though  the  tip  speed  is  slightly  less 
than  the  previous  case.  The  results  show  good  agreement  of  the  calculations  with  the  experimental 
data,  especially  near  the  tip. 

Figure  5  shows  the  extent  of  shock-induced  boundary  layer  separation  for  the  transonic  cases 
discussed  above.  These  are  delineated  as  surface  particle  flow  details  and  are  created  by  releasing 
fluid  particle  tracers  at  one  grid  point  above  the  surface  and  forcing  them  'o  stay  in  that  plane.  Such 
a  view  mimics  the  surface  oil  flow  details  measured  in  a  laboratory  experiment.  Figure  5a  shows 
the  details  of  this  flow  for  the  case  of  Mup  =  0.877  and  9  =  8°.  The  separation  and  reattachment 
locations  are  apparent  in  this  figure.  It  is  seen  that  this  flow  condition  produces  a  mild  shock- 
induced  separation  in  the  outboard  part  of  the  blade.  In  contrast,  the  shock-induced  separation 
and  viscous-inviscid  interaction  are  much  stronger  for  the  case  of  Mup  =  0.794  and  9  =  12°.  The 
surface  panicle  flow  pattern  for  this  more  severe  case  is  shown  in  figure  5b.  As  seen,  the  extent  of 
the  separation  is  much  larger  for  this  flow  condition  than  for  the  case  of  figure  5a.  It  is  interesting, 
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however,  to  now  that  the  separation  patterns  in  the  tip  region  arc  approximately  the  same  for  these 
cases. 

A  general  comparison  of  the  present  results  with  the  experimental  data  can  be  made  by  exam¬ 
ining  the  bound  circulation  distribution.  Figure  6  shows  such  a  plot  of  dimensionless  circulation, 
r/fl  R1 ,  as  function  of  r/R  for  8  *  8*  case  and  tip  speeds  of  0.44  and  0.877,  corresponding  to  the 
data  presented  in  figures  2  and  3. 1  lere  r  is  the  radial  distance  from  the  rotation  axis,  R  is  the  radius 
of  the  rotor,  O  is  the  constant  angulai  velocity  of  the  rotor,  and  f*  is  the  circulation.  The  integrated 
local  lift  values  are  used  from  both  the  coarse  and  fine  grid  calculations  to  compute  the  dimension¬ 
less  circulation  shown  in  figure  6.  Also  shown  are  the  integrated  data  from  (he  experiments,  which 
were  reported  to  be  essentially  independent  of  lip  speed.  The  calculations  show  a  fair  agreement 
with  the  experiments,  except  in  the  inboard  pan  of  the  blade.  This  suggests  that  only  the  near-field 
effects  of  the  (ip  vortex  arc  captured  as  well  as  desired.  There  are  two  possible  reasons  for  the 
poor  agreement  in  the  inboard  pan  of  the  blade.  First,  the  vortex  wake  becomes  diffused  in  the 
far- fie  Id  grid,  so  its  induced  effect  is  significantly  diminished.  Second,  the  inboard  plane  boundary 
condition  may  be  indcquatc.  In  contrast  to  the  experimental  observation,  the  present  calculations 
show  some  dependency  on  the  blade  tip  speed. 

In  the  tip  region  the  agreement  is  also  not  very  good.  This  may  be  due  to  the  bevel  lip  that 
is  used  in  the  computation  compared  to  the  square-tipped  blade  in  the  experiments.  Overall,  how¬ 
ever,  the  surface  pressure  distributions  appear  to  agree  better  with  the  experiments  than  the  bound- 
circulation  distribution.  Relatively  minor  discrepancies  in  the  pressure  distributions  near  the  lead¬ 
ing  edge,  where  the  experimental  transducer  locations  are  relatively  sparse,  seem  to  translate  into 
significant  differences  in  the  circulation  distribution. 

The  chief  advantage  of  the  Navicr-Stokcs  methods  is  to  predict  the  separated  flow  in  the  lip 
region  and  the  associated  detailed  structure  of  the  (ip  vonex.  The  prediction  of  the  overall  shed- 
wake  geometry  is  the  most  important  step  in  (he  process  of  accurate  modeling  of  the  complete  hover 
fiowfield.  The  ability  to  keep  this  shed  wake  (including  the  vortex  structure)  intact  from  diffusing 
due  to  the  numerical  dissipation  is  a  more  complex  issue.  The  ability  to  convect  this  shed  wake 
without  numerical  dissipation  determines  if  the  inflow  in  the  inboard  pans  of  the  blade  is  correct. 

Figure  7  shows  a  near-field  view  of  the  lip  vonex  panicle  path  trajectory  for  the  experimental 
conditions  of  Mhp  -  0.794  and  0=12°  corresponding  to  figure  4.  These  trajectories  are  generated 
by  releasing  panicles  of  fluid  in  the  vicinity  of  the  tip  of  the  blade  on  both  surfaces  and  allowing 
them  to  move  freely  in  time  and  space.  It  is  apparent  from  this  that  the  panicles  released  right  on 
the  tip  become  braided  and  stay  together  in  the  vicinity  of  the  core.  As  observed  before  (ref.  3 1 ),  the 
process  of  formation  of  die  tip  vortex  involves  braiding  of  fluid  panicles  from  both  upper  and  lower 
surface  of  the  blade.  As  the  process  of  braiding  of  fluid  panicles  from  upper  and  lower  surfaces 
continues,  the  tip  vonex  lifts  up  from  the  upper  surface  and  rolls  inboard  in  the  downstream  wake. 

After  identifying  the  fluid  panicles  in  the  vicinity  of  the  core  in  figure  7,  fewer  panicles  were 
released  on  the  tip  of  the  blade  in  the  proximity  of  the  quarter-chord  region  to  trace  out  a  trajectory 
of  the  tip  vonex  path.  Figure  8  show  two  views  of  this  trajectory.  The  computed  tip  vortex  trajectory 
in  space  for  a  single  blade  is  shown  in  figure  8a.  Figure  8b  shows  a  view  of  the  tip  vonex  looking 
from  the  top  which  highlights  the  contraction  of  the  wake.  The  contraction  of  the  wake  at  180° 
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and  360*  azimuthal  positions  is  approximately  12.8%  and  18.2%  of  the  radium,  respectively,  in 
agreement  with  the  experimental  observation  of  12.5%  and  17%  for  this  (low  condition. 

Fine  Grid  versus  Coarse  Grid  Results 

The  results  presented  in  the  prcceeding  sections  were  calculated  on  a  fine  grid  of  nearly  one 
million  points.  The  initial  test  calculations  were  made  primarily  on  a  coarse  Navier-Stokcs  grid  of 
109x36x31  size.  This  grid  was  generated  by  removing  every  other  point  from  the  fine  grid  in  all 
three  directions.  The  outer  dimensions  of  the  grid  and  the  grid  topology  are  thus  the  same  as  for 
the  pre"ious  fine  grid. 

Figure  9  shows  a  comparison  of  surface  pressure  distributions  for  the  fine  and  coarse  grids 
for  the  experimental  condition  of  Afi ,>  »  0.877,  8  *  8*  and  Re  *  3.93x  10*.  ft  is  surprising  to 
see  such  good  agreement  of  the  coarse  grid  results  overall  with  those  of  the  fine  grid  and  with 
the  experiments.  In  the  regions  where  the  shocks  are  very  strong,  there  are  slight  differences  as 
expected.  The  results  inboard  of  y/R  *  0.50  show  a  bigger  difference  as  seen  from  figure  6.  These 
quasi-steady  results  for  the  coarse  grid  took  about  one  hour  of  CPU  (central  processor  unit)  time 
on  the  Cray-2  supercomputer. 


Euler  versus  Navier-Stokcs  Results 

As  discussed  earlier,  there  have  been  several  attempts  to  capture  rotor  wakes  using  Euler  meth¬ 
ods  (refs.  14-17).  The  vortex  formation  in  the  tip  region  of  a  wing  or  a  helicopter  Made  is  a  result 
of  complex  three-dimensional  separated  flow,  and  it  is  not  clear  how  the  Euler  methods  arc  able  to 
mimic  viscosity  and  separation  to  produce  a  vortex  structure.  Nevertheless  Jik\$c  Euler  methods 
have  been  able  to  predict  the  pressure  distributions  and  spanwise  blade  losing  reasonably  well 
for  the  outer  part  of  the  blade.  Against  this  background,  a  limited  comparison  of  surface  pressures 
has  been  made  for  the  Euler  and  Navier-Stokcs  methods  calculated  on  the  same  fine  grid  of  about 
one  million  points.  It  may  be  noted  that  the  Euler  version  of  the  code  did  not  exhibit  any  stability 
problems  with  this  fine  Navier-Stokcs  grid. 

A  typical  comparison  of  the  Euler  results  with  the  Navier-Stokcs  results  is  presented  in  fig¬ 
ure  10  for  the  experimental  test  condition  of  A {tip  -  0.877, 8  =  8“,  and  Re  =  3.93  x  106.  Because  it 
neglects  viscous-inviscid  interaction,  the  Euler  method  overpredicts  the  shock  wave  strength  and 
position  for  y/R  0.80.  Otherwise,  the  Euler  results  are  in  good  agreement  with  the  Navier-Stokcs 
results,  which  show  mild  shock-induced  separation  for  this  flow  condition  (sec  fig.  5a).  The  over¬ 
all  agreement  of  surface  pressures  certainly  does  not  reflect  the  details  of  the  flow  near  the  blade 
surface,  especially  the  separation  pattern  and  vortex  wake  details  as  predicted  by  the  Navier-Stokcs 
method.  The  details  of  the  wake  structure  need  to  be  investigated  further. 


CONCLUSIONS 


The  lifting  hovering  rotor  flowfield  is  calculated  by  means  of  an  implicit,  completely  upwind, 
finite-difference  numerical  procedure  for  the  solution  of  »hin  layer  Navier-Stokcs  equal!  ns  using 
a  cylindrical  C-H  grid  topology  and  body  fixed  coosxiinates.  Th,  vortex  wake  and  its  induced 
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effects  art  captured  as  a  pan  cf  the  overall  numerical  solution  by  the  use  of  a  periodicity  condition, 
and  the  method  therefore  does  not  use  any  ad  hoc  wake  models.  The  present  numerical  results 
are  in  good  agreement  with  the  experimental  data,  and  they  represent  an  improvement  over  the 
previously  published  Navicr-Stokes  results  that  used  a  simple  wake  model.  Therefore,  the  method  is 
promising.  However,  several  important  issues  such  as  drag,  power,  and  the  detailed  wake  geometry 
remain  to  be  examined  in  detail. 

The  good  agreement  of  the  surface  pressures  predicted  by  the  Euler  method  with  those  of 
Navicr-Stokes  results  seems  to  suggest  that  the  details  of  surface  flow  including  separation  and 
tip  vortex  details  are  not  important  for  predicting  airloads.  This  needs  further  investigation.  The 
robustness  of  the  present  methodology  for  Euler  calculations  is  also  demonstrated.  Comparison 
of  coarse  and  fine  grid  results  indicate  that  the  farfield  wake  effects  are  not  as  well  captured  with 
coarse  grids.  The  numerical  method  is  fairly  efficient  and  runs  at  14S  M  FLO  PS  on  the  Cray*2 
supercomputer.  The  quasi-steady  Navier-Siokes  calculations  presented  here  for  coarse  and  fine 
grids  took  approximately  1  hour  and  IS  hours  of  CPU  time,  respectively,  on  this  machine. 
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Farfield  Boundary 


Fig.  1  Coarse  C-H  cylindrical  grid  topology  for  a  two-bladed  rotor;  a)  view  in  the  plane  of  the 
rotor,  and  b)  isometric  view  showing  the  grid  boundaries  for  a  single  blade. 


Fig.  2  Comparison  of  surface  pressures  for  a  lifting  hovering  rotor;  Mup  =  0.44,  9  =  8°,  and 
Re  =  1.92x10s. 
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-  Navler-Stokes  captured  wake  results  -  Present 

- Havior-Stokes  prescribed  wake  results  -  Re(.  9 

•  O  Experimental  data  -  Rcl.  29 


Fig.  3  Comparison  of  surface  pressures  for  a  lifting  hovering  rotor;  Mt,p  =  0.877,  0  =  8°,  and 
fle*  3.93x10s. 
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(a) 


Blade  Motion 


k  Leading  Edge 


(b) 


Fig.  5  Computed  surface  panicle  flow  detail  highlights  the  shock-induced  boundary  layer  separa¬ 
tion  for  the  flow  conditions  of  a)  A fj,p  =  0.877, 0  =  8°,  and  Re  =  3.93  x  10 6 ,  and  b)  A/<ip  =  0.794, 
0=  12°,and  J?c  =  3.55xl06. 
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Navler-Stokes  fine  grid, 
Navler-Stokes  coarse  grid 
Navler-Stokes  fine  grid 
Navler-Stokes  coarse. grid 


]  Mt|p  =  O.fi 
]Mtlp=0.44  J 


>  Present 


Pig.  6  Comparison  of  bound  circulation  distribution  for  the  case  of  collective  pilch  0  =  8°  with  tip 
speeds  of  Mt,p  =  0.44  and  0.877. 


Fig.  7  Calculated  tip  vortex  panicle  flow  details  showing  the  near- field  view  for  the  condition 
Mlip  =  0.794, 0  =  12°,  and  Re  =  3.55x  10  6 . 
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-  Navlcr-Stokes  captured  wake  results  -  Fine  grid  (217  *  71  *  61) 

- Navier-Stokcs  captured  wake  results  -  Coarse  grid  (109  *  36  *  31) 

•  O  Experimental  data  -  Rel,  29 


Fig.  9  Comparison  of  surface  pressures  with  coarse  and  fine  grids  for  the  case  of  Mtip  =  0. 
0  =  8°,  and  Re  =  3.93x  IQ6. 
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